Helium line formation and abundance in a solar active region 
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X ' ABSTRACT 

An observing campaign (SOHO JOP 139), coordinated between ground based and SOHO 
instruments, has been planned to obtain simultaneous spectroheliograms ol the same active region 
in several spectral lines. The chromospheric lines Ca II K, Ha and Na I D as well as He I 10830, 
5876, 584 and He II 304 A lines have been observed. The EUV radiation in the range A < 500 A 
and in the range 260 < A < 340 A has also been measured at the same time. These simultaneous 
observations allow us to build semi-empirical models of the chromosphere and low transition 
region of an active region, taking into account the estimated total number of photoionizing 
photons impinging on the considered active region and their spectral distribution. We obtained 
a model that matches very well all the observed line profiles, using a standard value for the He 
abundance ([He]=0.1) and a modified distribution of microturbulence. For this model we study 
the influence of the coronal radiation on the computed helium lines. We find that, even in an 
active region, the incident coronal radiation has a limited effect on the UV He lines, while it 
results of fundamental importance for the D3 and 10830 lines. Finally we build two more models 
assuming values of He abundance [He]= 0.07 and 1.5, only in the region where temperatures 
are larger than 1. x 10 4 K. This region, between the chromosphere and transition region, has 
been indicated as a good candidate for processes that might be responsible for strong variations 
of [He] (?). The set of our observables can still be well reproduced in both cases changing the 
atmospheric structure mainly in the low transition region. This implies that, to choose between 
different values of [He], it is necessary to constrain the transition region with different observables, 
independent on the He lines. ^ 



Subject headings: Sun:activity — Sun:abundances — Sun:Helium 



1. Introduction 

Helium is a very peculiar element in the sense 
that, despite its high abundance in the universe, 
its photospheric solar abundance is unknown. Un- 
fortunately it is undetectable in the visible photo- 
spheric spectrum and, until recently, the so-called 
solar photospheric abundance has been derived us- 
ing theoretical stellar evolution models. A com- 
monly accepted value for the ratio [He] = N# e /N# 
is 0.1 and represents the helium abundance of the 
nebula from which the solar system formed. How- 
ever, there are now strong indications that this 
value is too large: the inversion of helioseismic 
data by different authors leads to values of [He], 
in the convection zone, in the range 0.078 - 0.088 
(for a review, see ?). A reduction of He abun- 
dance of roughly 10% below its initial value could 
be explained including diffusion in the theoretical 
evolution models (??). 

The values quoted above refer to the outer con- 
vection zone and to the photosphere. Measure- 
ments of [He] in the solar wind and interplanetary 
medium exhibit a high variability, from almost 
negligible values (« 0.01) in equatorial regions, to 
high values (« 0.3) in transient events. An aver- 
age value in the fast component of the solar wind 
is about 0.05, thus about a factor of two lower 
than the photospheric [He] (??). This value is also 
in contrast with the general pattern of abundance 
variations in the fast wind, where elements with 
First Ionization Potential (FIP) > 10 eV typically 
retain their photospheric abundance, as opposed 
to low FIP elements (< 10 eV) which appear en- 
hanced (for a review, see ?). 

Direct measurements of [He] in the solar atmo- 
sphere closer to the surface are available only for 
the corona. In the quiet corona, at R ~ 1.1 Rq, 
? and ? found [He] w 0.08, suggesting no helium 
depletion in the low corona, while in a equatorial 
streamer, at R « 1.5 Rq, ? found an upper limit of 
0.05. More recently, ?? have estimated an abun- 
dance of 0.05 and 0.04 respectively, at distances 
R < 1.1 R Q . 

Apart from the observed variability and/or in- 
consistencies of the measurements in the corona 
and the solar wind, the general pattern of the FIP 
effect suggests that elemental fractionation should 
occur in the chromosphere, at T < 10 4 K (?). De- 
tailed theoretical calculations including diffusion 



effects (??) seem to confirm that, indeed, [He] un- 
dergoes strong changes from the chromosphere to 
the lower corona, before stabilizing at the values 
observed in the solar wind. 

Therefore it would be very important to de- 
termine the He abundance at chromospheric and 
transition region (TR) levels directly from line 
profiles. The strongest He lines in the solar at- 
mosphere are in the extreme ultraviolet (EUV) 
below 584 A. From the ground, the only observ- 
able He line in the quiet solar atmosphere (on 
the disc) is in the near infrared at 10830 A, but 
when active regions or flares are observed the He 
lines are enhanced and in particular the D 3 line 
at 5876 A can also be detected. The excitation 
of these subordinate lines is at least partially due 
to a photoionization-recombination mechanism or 
"PR" process (??), i.e. to direct photoionization 
of chromospheric helium atoms by coronal EUV 
radiation shortward of 504 A and successive re- 
combination to the excited levels of He I. Observa- 
tions actually suggest that this process might not 
be very important in quiescent regions, at least 
for the strong resonance, hydrogen-like He II lines 
such as the Lyman-a at 304 A (e.g.: ???). How- 
ever, no conclusive determination of the role of this 
process is available yet for active regions, where 
the coronal EUV radiation field is much more in- 
tense than in quiet areas. Hence, the study of the 
solar He spectrum requires not only spectral ob- 
servations of an active region (or flare) in a very 
large spectral range (from EUV to near infrared), 
but also an estimate of the coronal EUV radiation 
impinging on the observed active region. 

To this aim we planned an observing campaign 
(SOHO JOP 139) coordinated between ground 
based and SOHO instruments to obtain simultane- 
ous spectroheliograms of the same area in several 
spectral lines, including four He lines, that sam- 
ple the solar atmosphere from the chromosphere 
to the transition region. During this campaign we 
observed the region NOAA 9468 (cos#=0.99) on 
May 26, 2001 from 13:00 to 18:00 UT. A small 
two- ribbon flare (GOES class C-l), developed in 
this region around 16 UT and its dynamics during 
the impulsive phase has been already studied in ?. 

In this paper we concentrate on the pre-flare 
phase. We present semi-empirical models of the 
atmosphere of the active region (at 15:20 UT) con- 
structed to match the observed line profiles from 
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the chromosphere to the TR and taking into ac- 
count the EUV radiation in the range 1 - 500 A. 
The comparison between the observed and com- 
puted He line profiles allowed us to test the rele- 
vance of the PR process and of the He abundance 
values in an active region. In a following paper we 
will study the formation of He lines in the flare. 

2. Observations 

A detailed description of the observing program 
is given in ?. We recall here that spectrohcli- 
ograms were acquired with the Horizontal Spectro- 
graph at the Dunn Solar Telescope (DST) of the 
National Solar Observatory / Sacramento Peak in 
the chromospheric lines Ca n K, Ha and Na I D 
as well as in the He I lines at 5876 (D 3 ) and 10830 
A. The full field of view of the DST, 170" x 170", 
was covered in about 5 minutes, with a sampling 
step of 2". Correcting for offsets among the differ- 
ent detectors resulted in a final useful field of view 
(FOV) of 160" x 140" with an effective resolution 
of 2". 

During the same period, spectroheliograms of 
the active region were obtained in the spectral 
windows around the He I 584 A and the He II 
304 A spectral lines with the Normal Incidence 
Spectrometer (NIS) of the Coronal Diagnostic 
Spectrometer (CDS) aboard SOHO. The slit was 
stepped 6" covering a 148" wide area in ~ 5.5 min- 
utes. The final useful FOV was 148" x 138" with 
an effective resolution of 6" x 3.4". Ground-based 
and CDS data were aligned using SOHO/MDI im- 
ages as a reference. We estimate the error around 
few arcseconds. CDS and ground based spectra 
are simultaneous within 2 minutes. 

The EUV flux was monitored by the Solar 
EUV Monitor (SEM) instrument aboard SoHO 
(?). The CELIAS/SEM instrument provides cali- 
brated total photon counts in the range A < 500 A 
(zero-th order, SEM ), and in the range 260 < A < 
340 A (first order, SEMi), at 1 AU. 

Finally, the Extreme ultraviolet Imaging Tele- 
scope (EIT; ?) provided synoptic series of full disk 
images centered on 171, 195, 284 and 304 A around 
13:00 and 19:00 UT. 



3. Derivation of the photoionizing radia- 
tion 

In Sec. 1 we have mentioned the role of EUV 
radiation, emitted mainly by coronal plasma, in 
the excitation of helium lines. The relevant radi- 
ation we need to estimate falls in the wavelength 
range below the photoionization threshold of He I 
at 504 A (the photoionization threshold for He II 
is at 228 A). A misconception, commonly found 
in the literature, is that it would be sufficient to 
have an estimate of the intensities just near the 
photoionization thresholds, because of the sharp 
decrease with wavelength of the photoionization 
cross-sections - which are proportional to ss A 2 
and to ~ A 3 for He I and He II, respectively. How- 
ever, as discussed more in detail by ?, even pho- 
tons of considerably shorter wavelength can effi- 
ciently photoionize helium atoms and ions: what 
matters most is the ratio between helium and hy- 
drogen photoionization cross sections, at least for 
A > 20-50 A (at shorter wavelengths, inner-shell 
photoionization of metals absorb most of the pho- 
tons). The wavelength of the photoionizing pho- 
tons mainly determines how deeply in the atmo- 
sphere they can penetrate, since the photoioniza- 
tion absorption cross-sections decrease with wave- 
length. This point will be further illustrated in 
Sec. 4.3. 

Therefore, for the calculation of the He line 
profiles we needed to estimate the total number 
of photoionizing photons impinging on the tar- 
get active region, plus their spectral distribution. 
We obtained reliable estimates of these quanti- 
ties combining both observational constraints from 
CELIAS/SEM and EIT, and theoretical tools. 
The adopted procedures are illustrated in the fol- 
lowing sections. 

3.1. Coronal irradiance 

The solar irradiances measured at about 15:20 UT 
in the SEM and SEMi wavebands are $ SEM o = 
5.05 x 10 10 photons cm" 2 s" 1 and $ SEM i = 
2.57 x 10 10 photons cm~ 2 s _1 , respectively, with 
noise-like variations of the order of 0.1%. Even 
the subsequent Cl.l-class flare, starting at 16:01 
UT (?), produced variations in the SEM irradi- 
ances limited to at most 1%. Both quantities 
include the contribution of the He II 304 A line, 
while $ SEM o contains also the other He II reso- 
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nance features. Since the He lines and continua 
will be self-consistently computed in our radiative 
transfer calculations, we need to exclude those 
contributions to obtain the input values for our 
model calculations. 

The He II 256 A line has been observed to be at 
least 15 times weaker than the 304 A line in quiet 
regions (?). For an active region, observations of 
the first four terms of the resonance series indi- 
cate that the He II 304 A is at least 40 times more 
intense than the 256 A line, and about 25 times 
more intense than the sum of the 256, 243, and 
237 A terms (?). Evaluating the contribution of 
the continuum at A < 228 A is considerably more 
difficult (see, for instance ?). In an optically thin 
recombination spectrum, the contribution of that 
continuum could be as high as 30% of the total 
number of photons. However, the existing obser- 
vations of the resonance series (such as those cited 
above) strongly suggest a formation of the He II 
spectrum at high optical depths. Therefore, we 
may expect a considerably smaller contribution of 
the continuum and of other lines, as already sug- 
gested by ? and ?. Overall, we thus estimate that 
the contribution of the He II features below 256 A 
is less than 10% of that of the He II 304 A line, so 
that only the latter line need be taken into account 
in the measured SEM irradiances. 

The contribution of the He II 304 A line (plus 
the nearby Si XI 303 A line) to SEMi irradi- 
ances has been estimated by ? at approximately 
50%-60% for observations acquired during 2001. 
Adopting an average value of 55%, $ Hc + Sl = 
1.4 x 10 10 photons cm~ 2 s _1 (the contribution 
of the Si XI 303 A line is about 10%). Sub- 
tracting this value from the SEMo irradiances, we 
obtain the solar irradiance due to coronal lines 
only: $ COTOna « 3.64 x 10 10 photons cm~ 2 s" 1 . 
The uncertainty in this value is mainly due to 
the uncertainty in the absolute calibration of the 
He II 304 A line relative to the other lines in the 
CDS spectra used by ?. For CDS, the maximum 
uncertainty of the 2 n< ^ order calibration, relevant 
for the He II 304 A line, can be conservatively 
estimated at 30% (?). Thus, the above value 
of $ Hc + Sl has a maximum relative error of about 
14%, which translates into a maximum error of 
0.2 x 10 10 photons cm" 2 s" 1 for $ COTOna . If wc 
assimilate this maximum error to a 3-cr uncer- 
tainty, the standard deviation for <j> corona C an be 



estimated to be 0.06 x 10 10 photons cm 2 s 1 . 
We thus finally obtain: $ COTOna = (3.6±0.1) x 10 10 
photons cm~ 2 s _1 (1-cr error). In the estimated 
uncertainty for this value, we can safely include 
the errors introduced by neglecting the contribu- 
tions of the Si XI 303 A lines and of the oth- 
ers He II features, contributions which should be 
roughly of the same order, as mentioned above, 
but of opposite sign. 

3.2. Total number of ionizing photons at 
the solar surface 

To obtain the contribution from the area 
around the DST slit to the total EUV solar ir- 
radiance we used the spatial distribution of EUV 
emission provided by EIT images. The details of 
the method are described in the Appendix and 
we illustrate here only the main assumptions and 
results. Basically, to obtain the mean value of 
in any region of the EIT FOV, 
we scale the SEM calibrated irradiances, <l> SEMo 
and $ SEMl , by the ratio of the mean counts in the 
region over the total counts in the full EIT FOV 
(Eq. A4). 

This procedure assumes that the total counts 
over the detector in an EIT band (after basic 
CCD processing, including flat-fielding) are pro- 
portional to the irradiance, <I> EIT , in the same 
band. This is true if the contribution to $ EIT 
from the corona beyond 1.2-1.4 Rq (the limit of 
an EIT image) is negligible. In fact, inspection of 
the EIT images used in our analysis reveals that 
only a small fraction of counts comes from areas 
above 1.1 i? Q , and most of those counts are prob- 
ably due to light scattered in the telescope. 

The strongest assumption in this procedure 
however, is that the images in the EIT wavebands 
are good linear proxies of the intensity integrated 
within the SEM wavebands. In the case of the 
A < 500 A range (SEM ), all EIT bands (around 
171, 195, 284, and 304 A) include strong lines 
which contribute significantly to the total flux in 
that range. Thus, we may expect such an assump- 
tion to be quite reasonable. The case of the SEMi 
range is more straightforward: the EIT 304 wave- 
band, although narrower, is centered at the same 
wavelength (« 300 A) , and we may expect EIT 304 
images to be good proxies of the spatial variation 
on the disk of the radiance in the SEMi waveband. 



4 



On the other hand, the assumption of linear 
correlation of EIT intensities with integrated in- 
tensities in the SEM wavebands, is equivalent to 
assume that the spectral distribution of EUV pho- 
tons in those wavebands does not change with the 
location on the Sun and that only the overall in- 
tensity of the lines changes. This assumption is 
consistent with the calibration procedure of the 
SEM data (?), and is estimated to affect SEM 
measurements by about 10% (?). 

Hence, we considered a region of (70") 2 cen- 
tered at the position of the DST slit (solar rotation 
being taken into account), and we computed the 
conversion factors from <j> SEM to the mean value of 
jSem usm g ^ w0 gets of EIT synoptic images taken 
before and after the DST and CDS observations 
(at 13:00 and 19:00 UT). Since we already sub- 
tracted the contribution of the He II 304 A line 
to the SEM irradiance, we did not consider the 
EIT 304 images. In this way, we obtained an aver- 
age conversion factor c = 3.3±0.6 x 10 4 sr _1 , where 
the uncertainty was obtained from the standard 
deviation of the distribution of counts within the 
region of interest. Therefore, the intensity of coro- 
nal lines, integrated in the A < 500 A wavelength 
range and averaged over the area of the DST slit 

ig jcorona = ^corona x C = (1.2 ± 0.2) X 10 15 pllO- 

tons cm~ 2 s _1 sr _1 . The uncertainty is obtained 
combining the uncertainties of $ corona and c. 

In order to validate this approach, we compare 
the intensity measured directly from CDS in the 
considered area (/ Hc + Si = (1.1 ±0.1) x 10 15 ) with 
the one obtained from the SEMi irradiance us- 
ing the described procedure. The mean conver- 
sion factor from the two synoptic EIT 304 im- 
ages is (5.8 ± 2.5) x 10 4 sr" 1 and therefore / Hc + Si 
= $Ho+si x c = ( _ 8 ± 0.4) x 10 15 photons cm" 2 
s _1 sr _1 . The agreement between these estimates 
and the values measured with CDS is well within 
the uncertainties. 

3.3. Angular and spectral distribution of 
photoionizing photons 

The values computed in the previous section 
refer to intensities incident perpendicularly to the 
solar surface. In our following calculations, we as- 
sume that the photoionizing incident intensities 
are constant with direction. It is very difficult 
to evaluate this assumption, since the variation 
with angle of the photoionizing radiation depends 



on the spatial distribution of the coronal plasma 
around the region of interest. For instance, had we 
assumed a vertically stratified, plane-parallel, op- 
tically thin corona, we would have obtained an in- 
tensity varying as where [i is the cosine of the 
angle with the normal. Under such an assumption, 
the variation of the mean intensity, J, (and thus 
of the photoionization rate) with optical depth, 
r, is J(t) ~ Ei{t), where Ei(x) is the first ex- 
ponential integral. The constant intensity case we 
adopted gives instead a dependence J(r) ~ E 2 (t). 
However, at large depths the two assumptions give 
the same result: J(r) ~ exp(— t)/t. The differ- 
ences are strong only in the shallower layers, with 
the ratio E 2 {t)/E\ (t) falling below 2 already for 
r > 0.26. 

For the wavelength distribution of the ionizing 
photons, we adopt a synthetic spectrum obtained 
with the CHIANTI database (version 4: ?) with 
the Differential Emission Measure (DEM) of an 
"average" active region, from ?. In the calcu- 
lations, we used standard elemental abundances, 
with some values scaled to mimic the FIP effect 
(further details in the file sunJiybricLext . abund 
in the CHIANTI package), and a value P c /k — 
3 x 10 15 K cm~ 3 , where P c is the electron pressure 
and k the Boltzmann constant. 

We need to scale the synthetic spectrum (ex- 
cluding of course the He II lines and continua) by 
a factor 1.3 ± 0.3 to obtain the value of / corona 
found in Sec. 3.2. We note that such a linear scal- 
ing is equivalent to assuming that the shape of the 
DEM distribution of the plasma in our region of 
interest is the same as that of the distribution used 
to compute the synthetic spectrum. The resulting 
spectrum is shown in Fig. 1. 

4. Atmospheric models 

Our goal was the determination of semi- 
empirical models of the chromosphere and low 
transition region of the active region, providing a 
good match with the observed line profiles. The 
diagnostics we used for the chromosphcric part of 
the models were the profiles of the Ca n K, Ha, 
the Na I D lines, and the He I lines at 5876 A (D 3 ) 
and 10830 A, averaged over a region of 6" x 3.4", 
which corresponds to the effective resolution of 
CDS. For the transition region we used the CDS 
profiles of the He I line at 584 A and of the He II 
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Ly-a line at 304 A. 

4.1. General characteristics 

The modeling was done using the program 
PANDORA (?). Given a T vs. z distribution, 
we solved the non-LTE radiative transfer, and the 
statistical and hydrostatic equilibrium equations, 
and self-consistently computed non-LTE popula- 
tions for 10 levels of H, 29 of He I, 15 of Fe I, 9 of 
C I, 8 of Si I, Ca I and Na I, 6 of Al I and He II, 
and 7 of Mg I. In addition, we computed 6 levels 
of Mg II, and 5 of Ca II. 

More detail on the modeling, and on the differ- 
ent assumptions and their validity can be found 
in ?. Here, we just point out that the Ca II 
line-profiles are computed in Partial Redistribu- 
tion and that line-blanketing is included consid- 
ering Kurucz's opacities. The atomic models we 
used for H and Ca II are described in ? and ?. For 
He I, we used the 29 levels model described in ? 
which, apart for the number of levels, differs very 
slightly from the model of ? The different compo- 
nents in the lines were treated assuming that the 
sublevels are populated according to LTE, as ex- 
plained in ?. For He II we used the 6 levels model 
of?. 

The calculations include incident radiation 
from the coronal lines, which affects in particu- 
lar the ionization balance in the H, He I and He II 
continua. As explained in Sect. 3.3, this radiation 
was obtained scaling by a factor of 1.3 the dis- 
tribution given by the CHIANTI database for an 
active region. As a starting point for the helium 
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Fig. 1. — Estimated spectral distributions of pho- 
tons around the DST slit, derived as described in 
Sec. 3.3. The solid histogram shows the spectral 
distribution derived from a CHIANTI synthetic 
spectrum; the grey area indicates l-cr uncertain- 
ties. 



abundance, we used the standard photosphcric 
value [Hc]=0.1. The microturbulence followed the 
semiempirical distribution given by ? as a func- 
tion of total hydrogen density. 

As studied by Fontenla et al. (2002), the inclu- 
sion of ambipolar diffusion of helium can have an 
important effect in the construction of a model, 
changing the energy balance and therefore the T 
vs. z structure. However, it must be remarked 
that their models of the low transition region are 
theoretical. We have made some trial runs in- 
cluding ambipolar diffusion but keeping our atmo- 
spheric structure unchanged, and found that the 
emitted profiles remain unchanged. 

4.2. Model construction 

Given the different sensitivity of various spec- 
tral lines to modifications in different parts of the 
model, we can describe its construction as a step 
by step process. 

The source function for Ha is sensitive to the 
structure of the chromosphere from the tempera- 
ture minimum region to the base of the transition 
region, even if the line center is formed higher up 
in the atmosphere. The Ca n K line, in partic- 
ular its wings and the Ki minimum, are formed 
at the temperature minimum and in the low and 
mid-chromosphere. Therefore, we fixed the deeper 
part of the chromospheric structure using these 
two profiles, and cross-checked it with the profile 
of the Na D lines, which are much less sensitive 
to these atmospheric layers. None of these lines 
depend on the radiation field. 

The profiles of the He I lines at 5876 A (D 3 ) 
and 10830 A are formed in two distinct regions. 
Most of the radiation in both lines is originated in 
the photosphere, which in the quiet Sun results in 
a weak absorption line at 10830 A and no notice- 
able line at 5876 A. However, in an active region 
there is also a chromospheric contribution which 
depends mainly on the coronal EUV incident ra- 
diation, but also on the thermal structure of the 
high chromosphere/low transition region, between 
1. x 10 4 and 2.5 x 10 4 K. We used these lines, 
therefore, to fix this region as a second step in the 
construction of the model. 

Finally, both ultraviolet lines, the He I line at 
584 A and the He II line at 304 A, depend on 
the structure of the low and mid-transition region, 
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from 3. x 10 4 to 5. x 4 K for the 584 A line, and 
up to 1. x 10 5 K for the 304 A line. As men- 
tioned in the Introduction, their dependence on 
the radiation field is not yet assessed, and should 
therefore be investigated. Hence, we used these 
lines in a final step as diagnostics of the structure 
of the transition region. 

In Fig. 2 we show the distribution of the tem- 
perature as a function of column mass for the 
model that gives the best match to the observa- 
tions (dashed line). Also shown are the regions of 
formation of the different lines we used as diag- 
nostics for the model construction. 

In Fig. 3 we show the computed and observed 
line profiles for Ha, the Ca n K and Na D lines, 
which do not depend on the radiation field. All 
computed profiles were convolved with the instru- 
mental response. The bars in the Figure indicate 
the r.m.s. of the profiles averaged over the area 
observed by CDS. It can be seen that the agree- 
ment found is very good, well whithin the r.m.s. 

In Fig. 4 we compare the observed helium lines 
with the profiles obtained from this model, again 
convolved with the instrumental response. We 
want to stress that for the 584 A the instrumental 
effects are rather heavy, washing out completely 
the strong central self-absorption of the computed 
profile (see following Section and Fig. 6). The 
He II 304 A line profile (not shown) is very sim- 
ilar to the 584 A line. For all the helium lines, 
although the central intensities agree with the ob- 
servations within the variability between pixels, 
the computed profiles result too broad. Further- 
more, the observed profile of the 10830 A line 
distinctively shows the weak blue component of 
the triplet, which in the computed profile appears 
washed out. 

The main parameter that contributes to the 
broadening of these lines is the microturbulent ve- 
locity (vt). As described above, in this model we 
adopted the turbulence distribution given by ?. To 
obtain a better match with the observed He lines 
we then tried a different microturbulence distribu- 
tion, changing v t from 10 to 5 km s _1 in the region 
with temperature ranging from 1. x 10 4 to 21. x 10 4 
K. Since the microturbulent velocity is included in 
the hydrostatic equilibrium equations, changing it 
implies changing the densities, and therefore the 
atmospheric structure. Hence, the whole proce- 
dure was repeated, until a satisfactory match with 
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Fig. 2. — Temperature vs column mass distribu- 
tion of the computed atmospheric models. The 
models shown were obtained with the microturbu- 
lent velocity of ? (dashed line), and with a value 
of 5 km/s in the region of formation of the lines 
(solid line, see text). The plage model P (?) is 
displayed for reference (dotted line). Also shown 
are the heights of formation of the different lines 
we used to build the models. The He I 584 line 
is formed at temperatures between 2.4 x 10 4 and 
3.4 x 10 4 K, and the He n 304 line is formed be- 
tween 9.5 x 10 4 and 1. x 10 5 K. 
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Fig. 3. — Observed (dashed line) and computed 
(solid line) profiles for Ha and the Ca II K and 
Na D lines. The bars indicate the r.m.s. of the 
averaged profiles (see text) and point out the vari- 
ability between pixels. 



the observations was found, and a new model was 
built. This model is also shown in Fig. 2 as a 
solid line, and the corresponding profiles for the 
He lines are shown in Fig. 4. The new profiles for 
Ha, the Ca II K and Na D lines (not shown) are 
very similar to the ones obtained with the former 
model. It can be seen that the match to He pro- 
files is now strongly improved. We have therefore 
adopted this model as our standard model for the 
observed active region. 

Finally, we checked a posteriori that the value 
of P e /k = 3 x 10 15 K cm~ 3 , assumed to com- 
pute the spectral distribution of the coronal radia- 
tion with the CHIANTI database, agrees with that 
found with our models. The most external point 
of the model for which we have an observational 
constraint, has a temperature T= 1.1 x 10 5 K and 
an electron density N e = 2.3 x 10 10 cm -3 , resulting 
in a very good agreement with P e /k — 2.5 x 10 15 
K cm- 3 ' . 

4.3. Effect of the EUV radiation 

To assess the influence of the coronal radiation 
on the profiles of the helium lines, we performed 
new calculations for different intensities of the ra- 
diation field, and for different spectral distribu- 
tions. Since the most important change in the 
profiles with the radiation field is in the central 
depth of the D 3 and of the 10830 A line, in Fig. 5 
we show the intensity in the center of the 10830 
line as a function of the intensity / corona . 

Each curve in the figure represents a different 
spectral distribution. In each case, we considered 
a step function, which is non-zero only in a limited 
spectral range, and zero outside it. In this way, we 
tested distributions with all the photons localized 
between 25 and 50 A, between 50 and 100 A, 100 
and 200 A, 200 and 300 A, 300 and 400 A and 
between 400 and 504 A. 

Several important conclusions can be drawn 
from this Figure. First, photons at short wave- 
length (A < 50 A) have a small influence on the 
line profiles. Second, for the same value of J corona j 
photons between 50 and 100 A are the most effec- 
tive in increasing the line central intensity. Third, 
the details of the spectral distribution between 100 
and 504 A are of little relevance for the resulting 
central intensity. 

Such effects, anticipated by ? and reprised in 
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Fig. 4. — Dotted line: Observed profiles for three 
helium lines. Dashed: computed profiles (con- 
volved with instrumental response) for the model 
adopting the microturbulcnce of ?. Full line: same 
as dashed, for our standard model where microtur- 
bulence has been reduced (see text). In the upper 
and middle panel, bars indicate the r.m.s. of the 
averaged profiles and point out the variability be- 
tween pixels. For the shallow D3 line (maximum 
central depth « 5%) the variability is very large 
and the bars are clipped at the limits of the plot. 



Sect. 3, can be understood in terms of the de- 
pendence on wavelength of the fraction of EUV 
photons effectively capable of photoionizing he- 
lium atoms and ions. The first effect is essentially 
due to the fact that short-wavelength photons are 
almost entirely absorbed by metals, via inner-shell 
photoionizations. At longer wavelengths, the rele- 
vant parameter is the ratio of helium and hydrogen 
photoionization cross sections. The fact that the 
ratio between 50 and 100 A is about 3-5 times the 
same ratio at 500 A, is consistent with the second 
conclusion drawn from Fig. 5. Likewise, the third 
result can be understood by noting that the ra- 
tio between the He I and hydrogen cross-sections 
changes in a relatively slow fashion at long wave- 
lengths. 

In light of these results, we note that, in prin- 
ciple, any given central intensity of the 10830 line 
could be obtained with different values of / corona j 
but only if assuming very different spectral distri- 
butions. For example, the observed central inten- 
sity can be obtained with 8 x 10 14 photons cm~ 2 
s~l sr _1 between 50 and 100 A, or with almost 
twice as many photons between 300 and 504 A. 
However, for a "stationary" active region such as 
in our case, the hypothesis of a dominant contribu- 
tion at shorter wavelengths is rather unplausible, 
so that our assumption of a synthetic spectrum 
obtained with the CHIANTI database (Sect. 3.3) 
seems fully justified. 

In Fig. 6 we show the effect of the coronal radia- 
tion on the 584 A profile. We compare the profiles 
for the 584 A obtained with no incident coronal ra- 
diation with those obtained with / COTOna = 1 x 10 15 
photons cm~ 2 s~l sr _1 with two different spectral 
distribution. We see that the profiles are very sim- 
ilar and that the maximum difference (in the peak 
intensity) is about 10%. The profile variation de- 
pends more on the spectral distribution than on 
the intensity of coronal radiation. As expected, 
when all photons are between 50 and 100 A and 
penetrate deeper in the atmosphere, the line does 
not change appreciably. A small change occurs in- 
stead when all photons are between 400 and 504 
A. We obtained a similar result for the He 11 304 
A line. 

Our analysis hence shows that, similarly to the 
quiescent case, also for active regions the coronal 
radiation has a much stronger effect on the chro- 
mospheric He lines than on the TR He lines. 
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4.4. Helium Abundance 
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Fig. 5. — Central intensity for the 10830 A line as 
a function of the intensity / corona for different step 
distributions of the field. 
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Fig. 6. — 584 A profile computed with no incident 
coronal radiation (full line), compared with the 
profiles computed for / corona = l xlO 15 photons 
cm~ 2 s _1 sr -1 but with different spectral distri- 
bution. Dashed line: all photons between 400 and 
504 A. Dotted line: all photons between 50 and 
100 A. 



In the previous Sections, we have shown how 
a semi-empirical model employing the standard 
value of [He] =0.1 can reproduce very well the ob- 
servational data. On the other hand, the actual 
value of [He] in different parts of the solar atmo- 
sphere is a matter of discussion (see Introduction), 
so it might be of interest to investigate how differ- 
ent values of [He] affect our computations. 

Changing the abundance modifies the total den- 
sity for a given height, and therefore affects the 
hydrostatic equilibrium. Hence, modifying [He] in 
the whole solar atmosphere would result in a com- 
pletely different atmospheric structure. In partic- 
ular, the resulting photosphere could not be con- 
strained by our observations. Furthermore, the 
region between the chromosphere and transition 
region has been indicated as a good candidate for 
processes that might be responsible for strong vari- 
ations of [He] (?). For these reasons, we changed 
[He] only from the point where T = 1. x 10 4 K 
outwards, where the He lines are formed, and re- 
computed the models until we found a satisfactory 
match with the observations. We built two differ- 
ent models, assuming values of the abundance of 
[He]=0.07 and [He]=0.15. 

In Fig. 7 we show these models, together with 
the one obtained for the standard abundance 
[He]=0.1, in the region above log(m)=-3, where 
they differ. In Fig. 8 we compare the 10830 and 
584 A profiles obtained with these three models. 
In general, the differences in the profiles arc very 
small, although a larger value of abundance corre- 
sponds to narrower profiles in both lines. This is 
due to the fact that to compensate for a smaller 
abundance the region of formation of the lines 
needs to be extended, to obtain enough helium 
atoms. This larger region between 1. x 10 4 and 
2.5 x 10 4 K can be noted in Fig. 7, and results 
in a broader line. Although the ratio between 
the peaks and the central absorption of the 584 
A changes by a factor of two with the extreme 
abundance values, the convolution with the broad 
CDS response destroys any difference. Only ob- 
servations of EUV lines at much higher spectral 
resolution might then help in discriminate between 
different abundance values. 
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5. Discussion and conclusions 

We obtained a large set of simultaneous and 
cospatial observations of an active region, includ- 
ing the chromospheric lines Ca n K, Ha and Na I 
D as well as the He I lines at 5876 (D 3 ), 10830 and 
584 A and the He II line at 304 A. The EUV radi- 
ation in the range A < 500 A and in the range 
260 < A < 340 A has also been measured at 
the same time. Using the program PANDORA 
we were able to build semiempirical atmospheric 
models to match this set of observables. The 
He lines and continua are self-consistently com- 
puted in the radiative transfer calculations. In 
particular we were interested in the role of the 
photoionization-rccombination process in the for- 
mation of the He lines and in the effect of possible 
variations of [He] in an active region. 

One of the main problems for the calculation 
of the He line profiles was to estimate the total 
number of ionizing photons impinging on the tar- 
get active region, and their spectral distribution. 
The adopted procedure combines the observations 
from CELIAS/SEM and EIT instruments for our 
region with the synthetic spectrum obtained with 
the CHIANTI database for an "average" active re- 
gion and a P c /k = 3x 10 15 K cm -3 . We needed to 
scale the synthetic spectrum by a factor 1.3 ± 0.3 
to obtain the estimated value of J COTOna (Sec. 3.2). 

We have two independent checks to validate 
this approach. First, we compared the He II 
304 intensity, measured directly from CDS in the 
considered area, with the one obtained from the 
SEMi irradiancc in the range 260 < A < 340 A 
and the agreement resulted well within the un- 
certainties. Second, a posteriori, we found that 
the electron pressure for the most external point 
of the model for which we have an observational 
constraint, agrees with the one assumed to com- 
pute the spectral distribution with the CHIANTI 
database. The agreement guarantees that the ob- 
tained model is self-consistent with the EUV coro- 
nal intensity used. 

A semiempirical model with a standard value 
for [He] and a modified distribution of microtur- 
bulence v t (reduced in the region of formation of 
the He lines) provides a good agreement with all 
the observations. However, we must stress that 
the broad CDS response function heavily conceals 
the details of the EUV line profiles, thus limiting 



their diagnostic relevance. 

For this model, defined as our standard model, 
we study the influence of the coronal radiation on 
the computed helium lines. We find that, simi- 
larly to the quiescent case, also in an active region 
the incident coronal radiation has a limited effect 
on the UV He lines, while it results fundamen- 
tal for the D 3 and 10830 lines. For these latter 
two lines, the photons of wavelengths below 100 A 
result very effective in increasing the line depth, 
confirming that for the calculation of the He line 
profiles it is crucial to correctly estimate the to- 
tal number of ionizing photons and their spectral 
distribution. While the assumption of an average 
distribution might not be critical in the case of a 
"stationary" active region, this could be not true 
in the case of a flare, where the coronal EUV ra- 
diation is probably harder. 

Finally, we tested how the helium abundance 
influenced our computed profiles, changing the 
value of [He] only in the region where the tempera- 
tures are larger than 1. x 10 4 K. We built two more 
models for [He] =0.15 and [He] =0.07 and found 
that also with these models we could match the 
observations as well as with the standard model. 
The differences in the computed lines are mostly 
evident for 584 and 304 A whose ratio between 
the peaks and the central absorption changes by a 
factor of two with the extreme [He] values. How- 
ever, given the coarse spectral resolution of CDS, 
they are not appreciable in our observations that, 
hence, do not provide enough constraints to choose 
between [He] values. Observations of the 584 and 
304 A lines with a spectral element smaller than 
0.025 AA/pix around 584 A might help in this 
issue. This resolution is currently achieved only 
by the SUMER spectrograph aboard SOHO. How- 
ever, SUMER has been restricted to off limb ob- 
servations since year 2000. Moreover it does not 
observe the He II Lyman alpha around 304 A. 
Thus we have to wait for a next generation EUV 
spectrographs aboard one of the fortcoming mis- 
sions such as Solar Orbiter. Furthermore, since 
the changes of the atmospheric structure with [He] 
are larger in the TR, it would be useful to con- 
strain this part of the atmosphere with other ob- 
servables, independent on the He lines. In par- 
ticular we suggest that simultaneous observations 
of EUV lines of C II, C in, Si II and Si in might 
be very important to define the structure of the 
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TR and hence help to discriminate between [He] 
values. 
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Fig. 7. — Temperature vs column mass distribu- 
tion of the atmospheric models obtained for dif- 
ferent values of the He abundance, in the region 
where they differ. Solid line: [He]=0.1, our stan- 
dard model. Dotted line: [He]=0.07. Dashed line: 
[He]=0.15. 
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Fig. 8. — Profile for two helium lines, for dif- 
ferent values of the helium abundance. Solid 
line: [He] =0.1, our standard model. Dotted line: 
[He]=0.07. Dashed line: [He]=0.15. For the 584 
A also the profiles before the convolution with the 
instrumental response are shown. 



A. Derivation of intensities from irradiance measurements 

If the Sun were a featurless disk, with no limb brightening nor off-limb emission, the disk-integrated flux 
(irradiance), <&, measured at a distance D would be easily related to the specific intensity, I, at the solar 
surface: 

$ = I x nR 2 /D 2 = Jxfl , 

where R and fl are, respectively, the solar radius, and the angular dimension of the solar disk (in stera- 
dians). 

In the case of a spatially resolved, full-disk solar image where it is possible to identify N discrete structures, 
each of area Si (i = 1 . . . N), and emitting a surface intensity 7j along the line of sight, the observed irradiance 
is: 

$ = S h x ^ = S h x 0i ' 

i i 

where fij = Si/D 2 is the angular size of the i-th structure. 

^From the above relation, the contribution to the irradiance from the individual structure i is therefore: 

Hence, the surface specific intensity in that structure can be written as 

/, = ^ = c,x$, (Al) 

where the conversion factor c, (in units of an inverse solid angle) is 

h 



(A2) 



The relations Al and A2 are useful if the total irradiance is available in absolute units (e. g. photons cm~ 2 
s _1 ), while a solar image in the same bandpass provides spatially resolved, uncalibrated (e. g. counts pixel -1 
s _1 ), intensities. It is also worth noting that the structures in Equations Al and A2 can be arbitrary: the 
intensities U can in fact be the intensities of a single pixel, or of a set of pixels. 

In the limit of continuous intensity distributions, it is easy to generalize those relations to obtain the mean 
intensity in an arbitrary region 1Z as 

"K J n 

Hence, the generalization of relations (Al) and (A2) is trivial: 

^calib. = $ calib. x (c)k; (A3) 

where the mean conversion factor is 

1 J K dO/ counts 

^ U = Shi ' Jdfl I counts ' (A4) 
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